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Abstract 

We perform a global analysis of neutron KJ 1 structure function data, obtained by 
combining proton and deuteron measurements over a large range of kinematics. From 
these data the lowest moments (n < 10) of the leading twist neutron structure func- 
tion are extracted. Particular attention is paid to nuclear effects in the deuteron, which 
become increasingly important for the higher moments. Our results for the nonsinglet, 
isovector p — n combination of the leading twist moments are used to test recent lattice 
simulations. We also determine the lowest few moments of the higher twist contribu- 
tions, and find these to be approximately isospin independent, suggesting the possible 
dominance of ud correlations over uu and dd in the nucleon. 



1 Introduction 



Inclusive lepton scattering experiments have for some time been the standard tool with 
which to study the internal quark-gluon structure of nucleons and nuclei. In unpolarized 
scattering, this structure is represented through two structure functions, F% and F 2 , 
which parameterize the differential cross section, 
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where OMott = (4a 2 E' 2 / 'Q 4 ) cos 2 | is the Mott cross section for the scattering from point 
particles, with Q the scattered lepton solid angle (9 is the polar scattering angle), and E' 
the recoil lepton energy. Here v is the energy transfer and M is the nucleon mass. The 
structure functions i^x 2 are functions of the four- momentum transfer squared q 2 = —Q 2 , 
and the Bjorken scaling variable x = Q 2 /2Mu. In the Bjorken limit (Q 2 ,v — ► 00, 
x fixed), the structure functions can be related to parton distribution functions, which 
give the light-cone momentum distributions of quarks (and gluons, or partons in general) 
in the hadron. At leading order, one has F 2 = 2xFi = x^ li e 2 {qi(x,Q 2 ) +q~i(x,Q 2 )). 
At finite Q 2 radiative corrections and higher twist effects give rise to a Q 2 dependence 
in the structure functions, leading to scaling violations. 

The connection between the partonic structure of hadrons and QCD can be most 
readily made through moments, or x-weighted integrals, of structure functions, which 
through the operator product expansion (OPE) can be related to hadronic matrix ele- 
ments of local operators. In the case of the F 2 structure function, to which we restrict 
ourselves in the following, the moment is defined by 

M n (Q 2 ) = I' dx x n - 2 F 2 (x,Q 2 ) . (2) 
J 

While in the Bjorken limit the moments are given by real numbers, independent of Q 2 , 
at finite Q 2 the renormalization group equations lead to a specific Q 2 dependence of the 
moments that can be calculated from perturbative QCD. In fact, observation of this Q 2 
dependence was instrumental in establishing QCD as the theory of the strong nuclear 
interactions. The moments can also be calculated nonperturbatively using lattice QCD 
techniques. 

A wealth of information has been accumulated over several decades on the structure 
of the proton, deuteron, as well as heavier nuclei. To completely describe the structure 
of the nucleon, however, requires knowledge of both the proton and neutron structure 
functions. Typically, the neutron structure function is obtained by combining data on 
the deuteron and proton. In most analyses the deuteron structure function is assumed to 
be given by a simple sum of free proton and neutron structure functions, F 2 —* F 2 +F 2 , 
which is a good approximation for some regions of kinematics (0.1 < x < 0.4). However, 
it is well known that at small x (x < 0.1) nuclear shadowing reduces the deuteron 
structure function relative to the free nucleon ^ 03 El El- while at large x (x > 0.4) 
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nuclear binding and Fermi motion corrections are expected to give rise to a small "nuclear 
EMC effect" in the deuteron. 

In terms of moments, any nuclear effect at large x will manifest itself in high moments 
(large n). It is crucial, therefore, that a moment analysis of the neutron structure 
function includes a careful treatment of nuclear corrections. In this paper we for the 
first time perform such an analysis, utilizing the recent high-precision data from Jefferson 
Lab together with earlier data on F^ an d F% ( see Refs. jSJE] and references therein). 

The earlier world data on the F2 structure function cover a wide range of kinematics, 
from Q 2 ~ to ~ 100 (GeV/c) 2 . However, the large- x domain has been almost com- 
pletely unexplored, due to either experimental limitations or intepretation difficulties. 
New, high precision data from Jefferson Lab have allowed one to fill this gap which has 
in the past prevented a complete moment based analysis of data. Jefferson Lab data 
cover the entire resonance region starting from the elastic peak (quasi-elastic for the 
deuteron) and extending to invariant masses of the produced hadronic system of about 
3.5 (GeV/c). This region is of particular importance for higher moments. In fact, al- 
ready for the n = 4 moment the contribution of this kinematic region is significant for 
Q 2 < 5 (GeV/c) 2 . As a consequence for Q 2 < 5 (GeV/c) 2 and n > 2 the uncertainties 
come mostly from Jefferson Lab data. 

The lowest moment (n = 2) is also sensitive to the small- x region, where much data, 
in particular for Q 2 < 5 (GeV/c) 2 , have been accumulated in previous experiments. For 
the proton F2 structure function the lowest x were achieved with the HERA collider 
data which reached down to x ~ 10~ 5 . For the deuteron F2 structure function, the 
minimum x values are somewhat larger, x ~ 10 -4 , and were obtained with muon beams 
in fixed target experiments at CERN and Fermilab. The extrapolation to x = has 
been carefully studied in Refs. Ej and the corresponding systematic uncertainties 
have been estimated. At Q 2 values above those currently accessible at Jefferson Lab 
{Q 2 > 5 (GeV/c) 2 ), the large-x extrapolation was also studied by a comparison of 
various models and the corresponding systematic uncertainties were included in the 
moment analyses. 

Unlike previous studies, which have focussed on the x dependence, and have ignored 
or underestimated the nuclear corrections, we extract the lowest few moments of the neu- 
tron F2 taking into account the systematic uncertainties associated with our knowledge 
of nuclear effects in the deuteron. Furthermore, this analysis directly uses the leading 
twist parts of the moments extracted from the F% and F^ data, which allows a more 
straightforward implementation of the impulse approximation. 

The general framework for analyzing the data is presented in Sec. [21 where we dis- 
cuss various models of nucleon light-cone momentum distributions in the deuteron. The 
extracted leading twist neutron moments are presented in Sec. |3 where we also discuss 
corrections to the simple one-dimensional convolution formalism, in the form of nucleon 
off-shell effects, as well as nuclear shadowing and meson exchange currents. The results 
for the extracted moments are compared with those obtained from several global parton 
distribution functions. Moreover we construct the nonsinglet, isovector p — n combina- 
tion of the leading twist moments and compare our results with available lattice ones. 
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As an interesting application of our extraction, in Sec. 0] we identify the isospin depen- 
dence of the higher twist contributions to the moments, and speculate about possible 
interpretations of our findings within a diquark picture. Finally, in Sec. we summarize 
our results and briefly outline some future directions. 



2 Neutron structure and nuclear effects in the deuteron 

In the absence of free neutron targets, information on the structure function of the neu- 
tron is typically extracted from measurements using deuterons, together with knowledge 
of the corresponding proton structure function. Since the deuteron is a nucleus, its nu- 
cleon constituents are bound, albeit weakly, and in general the structure of the bound 
nucleons will differ from the free nucleon structure. In this section we describe the var- 
ious nuclear corrections which enter in the extraction of the free neutron F2 structure 
function, and in particular its moments, from the structure function of the deuteron. 

We begin with the standard nuclear impulse approximation, in which the virtual 
photon scatters incoherently from individual nucleons in the deuteron, and discuss sev- 
eral implementations of the convolution formula. Following this we describe corrections 
to convolution, associated with nucleon off-shell effects, nuclear shadowing, and meson 
exchange currents. 



2.1 Impulse approximation 

For moderate to large values of x (x > 0.1) inclusive scattering off a deuteron can be 
described, using the optical theorem, in terms of the amplitudes j*N* — > "f*N* and 
N*d — > N*d, where N* refers to the fact that the nucleon is off its mass or energy shell. 
In this approximation interactions between the intermediate states of the j*N* and N*d 
subprocesses are neglected. Formally, the deuteron F2 structure function can then be 
written as 

F 2 D (x,Q 2 ) = jd 4 pTr[w(p,q)-A(p,p D )], (3) 

where p, pd and q are the virtual nucleon, deuteron and virtual photon four- momenta, 
respectively, the Dirac matrix W(p, q) = I Wo + f Wi + fj W2 is the truncated (off- 
shell) nucleon structure function, and A(p,po) = IAo+"f a Ai is the relativistic deuteron 
spectral function. The amplitudes Ao and A" depend on the deuteron wave functions. 

For an on-shell nucleon, the F% structure function is proportional to M Wo+M 2 Wi + 
p-q W2. In general the presence of several Dirac structures and corresponding "structure 
functions" Wj in the trace in the integrand in Eq. complicates the relation between 
the deuteron and nucleon structure functions. There are various ways which the formal 
expression in Eq. (j3J) can be simplified, however, as we now discuss. 



2.1.1 Nonrelativistic expansion 

The formal expression for the deuteron structure function F 2 D in Eq. (j3J) is Lorentz- 
invariant and can be evaluated in any reference frame. Without loss of generality, how- 
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ever, one can choose to work in the deuteron rest frame, and take the z-axis such that 
q z = — |q|. Performing a nonrelativistic reduction of the Dirac traces in Eq. Q, one 
can show that to order p 2 /M 2 in the Bjorken limit the product of the Dirac structures 
reduces to a single structure proportional to a generalized (off-shell) nucleon structure 
function [I], 

*?(*.« , ) = /^(i + S)l»i,(p)l , ^(f.fl'y) • w 

Here F 2 = ^(i^ + i 7 ^ 1 ) is the isoscalar nucleon structure function, z = (po+p z )/M is the 
light-cone momentum fraction of the deuteron carried by the interacting nucleon, where 
Po = M + ££> — p 2 /(2M) is the energy of the struck nucleon, and Ed = —2.22 MeV 
is the deuteron binding energy. The deuteron wave function contains the usual 
nonrelativistic S- and D-wave components, ipo and fa, 

^L|* D (P)| 2 = </>o 2 (IpI)+vI(IpI), (5) 

normalized such that 

— % y o d|p|p 2 i* D ( P )i 2 = i. (6) 

Note that in addition to the Bjorken variable and Q 2 , the nucleon structure function in 
® also depends on the nucleon off-shellness p 2 = p\ — p 2 . If one neglects the dependence 
of i<2 on P 2 , then Eq. (JIJ leads to the simple one-dimensional convolution formula in 
terms of the free nucleon structure function [3 IE] , 

N fM D /M 

F 2 D ^ v \x, Q 2 ) = / dz f°(z) F» (~,Q 2 ) , (7) 

where is the deuteron mass, and f®{ z ) ls the nonrelativistic nucleon (light-cone) 
momentum distribution in the deuteron, 

Physically the function f®(z) describes the Fermi motion and binding of nucleons in the 
deuteron. Note that some implementations of the convolution formula Q do not include 
the "flux factor" (1 +p z /M) jUJ. However, since this appears at C(p/M), consistency 
to order p 2 /M 2 requires its inclusion 4 . The wave function normalization © ensures 
that 

M D /M 

dzfg(z) = 1, (9) 

which guarantees that (for the valence component) the convolution term alone preserves 
the baryon number of the deuteron. 
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2.1.2 Relativistic effects 

Because the deuteron is a weakly bound nucleus, in most calculations of deuteron struc- 
ture the nonrelativistic approximation is adequate. In the large- x region, however, one 
is more sensitive to the large-momentum components of the deuteron wave function. In 
this case relativistic effects are expected to become more important. 

Since the distribution function in Eq. Q was obtained in the nonrelativistic 
approximation, it applies in the region \z — 1| <C 1. Outside of this region relativistic 
corrections may be important. In order to test the sensitivity to these effects we also 
consider a relativistic version of the distribution in Eq. (JHJ) which is calculated with 
the same nonrelativistic wave function but using relativistic kinematics, namely po = 
M D - E p , where E p = \J M 2 + p 2 is the energy of the spectator nucleon. Futhermore, 
one also replaces the ratio p z /M in the flux factor with p z /E p . 

While most deuteron wave functions calculated to date have been nonrelativistic, 
there have been some models which have treated the two-nucleon system relativistically 
jlUl 111! IT2| . In order to consistently utilize the relativistic wave functions, one needs 
to work with the full expression for the structure function in Eq. © directly. This 
necessarily involves modeling the off-shell nucleon structure functions, which introduces 
additional model dependence in the calculation [T3] . 

One can show, however, that the exact expression in Eq. © can be reduced to an 
on-shell convolution part, similar to that in Eq. (J7J), plus nuclear off-shell and relativistic 
corrections [TH H5]. 

F 2 D (x,Q 2 ) = F 2 D(conv \x,Q 2 )+5^F 2 D (x,Q 2 ) . (10) 

The convolution component here now involves a relativistic nucleon light-cone distribu- 
tion function |14j . 

-oo 16?r PO 

where E p = \J M 2 + p 2 = (M 2 - p 2 + M|>) /2M D is the energy of the (on-shell) nucleon 
spectator, and the upper limit of the nucleon virtuality p 2 in the integration is given by 
Pmax = zMp/2 — zM 2 /(2 — z). In addition to the usual S- and D-state wave functions 
u and w, the wave function ^/£>(p) also includes the relativistic P-states |1 1 1 112) 

^ |§ D (p)| 2 = ^o 2 (IpI) + ^l(IPl) + 1&(|P|) + ^L(IPl) , (12) 

where ipi s and ipu correspond to the singlet and triplet P-states, respectively. Note that 
the function /5 also includes the relativistic flux factor, (E p /po) z and therefore 
satisfies the same normalization condition as in Eq. ©, / dzf^ = 1. One should point 
out, however, that the separation of F 2 D into the two terms in Eq. 1)10(1 is not unique 
(the sum, of course, is unique). However, the choice (jll|) is the most natural one, and 
corresponds to that which has been followed in many phenomenological treatments. 
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2.1.3 Light-cone formulation 

An alternative to the covariant (instant-form) formulations of scattering from the deuteron 
is the Hamiltonian light-cone (LC) approach. Here the intermediate state nucleons are 
on-mass-shell, but off-energy-shell, so that one avoids the explicit p 2 dependence in the 
bound nucleon structure function in Eq. (jlj). 

Because the nucleons are on-mass-shell, the derivation of the one-dimensional con- 
volution from Eq. (|3J) is trivial since each of the three terms in the trace factor become 
proportional ^3]- The LC nucleon momentum distribution in the deuteron is given by 



where |p| m in = M\z — z(2 — z), and Vlj^p is the LC deuteron wave function. Since 
most of the two-nucleon phenomenology has traditionally been formulated within the 
instant-form approach, in practice for vj/^p one has usually employed the S- and -D-state 
wave functions as calculated in nonrelativistic potential models [see Eq. (jSjl]. 

Note that in the limit where the nucleons are free (i.e. not bound), the deuteron can 
be described in terms of an NN Fock state component. Inclusion of binding requires the 
addition of iViV-l-meson Fock states. Technically, binding or off-shell effects in the LC 
approach can be included by constructing a wave function which explicitly depends on 
an additional light-like four- vector related to the orientation of the null-plane JZJ • This 
introduces a number of additional wave functions, which would need to be included in 
/lc- Recently calculations of deuteron LC wave functions have been performed within 
this framework, and some of the additional wave functions found to be important at 
large momentum. While it is interesting to explore the phenomenology using these wave 
functions, in the present study we shall compare the results of the LC distribution ff c 
using the nonrelativistic S- and -D-state wave functions as used in the literature. 

In the next section we compare the various light-cone momentum distributions dis- 
cussed above, for different model deuteron wave function models. 

2.2 Nucleon momentum distributions 

Within the convolution approximation there are two sources of nuclear model depen- 
dence in the calculation of the deuteron structure function. Firstly, as discussed above, 
different approaches are used to formulate scattering from the two-nucleon system - 
nonrelativistic, relativistic, or light-cone. In addition, within each approach there is 
dependence on the deuteron wave function or AW potential. It is beyond the scope 
of this work to critically analyze the merits of the various models; instead we survey 
a reasonably large sample of models and estimate the systematic uncertainty from the 
spread in their predictions. 

In practice the deuteron wave functions are determined by fitting NN scattering 
data. Since most of the data are at low momentum, the wave functions are relatively 
well constrained at small p, but the large-p tail of the wave function is more model 
dependent. In terms of the light-cone variable z, the large-p region corresponds to 




(13) 
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Figure 1: Nucleon light-cone momentum distribution f D weighted by y 7 , where y — z/2, 
for various models: (a) nonrelativistic distribution (solid), distribution with relativistic 
kinematics (dashed), relativistic model 111)) (dotted), light-cone distribution all using the 

Paris wave {unction fiflf : (b) relativistic distribution ]11\) using the Bonn \lty (solid), Paris USf 
(dashed), and relativistic Gross wave functions from Refs. \11\[ (dotted) and \l^j (dot-dashed). 



large values of z. To emphasize the large- z region of f D we plot in Fig. ^the various 
nucleon momentum distributions, for different wave functions, weighted by y 7 , where 
y = z/2. Note that the momentum integration in each of the curves has been limited 
to |p| = 0.5 GeV (see below). As expected the model dependence is relatively weak for 
moderate y, and only becomes significant for y > 0.6. 

Uncertainties due to the nuclear models can also be quantified in terms of moments. 
This is actually more relevant for the present analysis since ultimately we extract the 
moments of the neutron structure function. If we define the n-th moment of the 
distribution f D {z) by 



then clearly the uncertainties in the large- z part of the distribution will be reflected in 
the uncertainties of the moments of f D at large-n, which are illustrated in Fig. [2j 

As mentioned above, the differences between the momentum distributions at large z 
arise because the deuteron wave function is not well constrained at large momentum. In 
fact, AW scattering data which have been used to constrain the deuteron wave function 
typically extend only to p ~ 0.4 GeV. At large momentum (p > 0.3 — 0.4 GeV) the 
inter-nucleon distances probed becomes smaller than the nucleon radius, so that here 
a description in terms of nucleon degrees of freedom alone may be questionable. To 
estimate the uncertainty arising from the large-p tail of the deuteron wave function, we 
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Figure 2: Moments (n = 2 — 10) of the nucleon light-cone momentum distribution f D (z) 
in Eq. i (a) for different deuteron wave functions: BonnCD 1201/ (solid), Nijmegen 93 \21f 
(dashed), Argonne v%s '221 (dotted), Gross 11 If (dot-dashed); (b) for the Gross wave function 
\12j with the momentum integration restricted to p max = 0.5 GeV (solid) and 1 GeV (dashed). 



examine the effect of truncating the momentum integration at p = p max - The moments 
of the distribution f D (z) with p max = 0.5 and 1 GeV are plotted in Fig. Hfb), where 
the distributions have been renormalized to satisfy Eq. @. For low moments (n = 2,4) 
the effect is negligible, while for the ti — 10 moment the effect of the Pm&x — 

0.5 (1) GeV 

cuts is 6% (1%). We estimate that the uncertainty introduced by our poor knowledge of 
the large-p components of the deuteron wave function is of the same order of magnitude 
as the wave function model dependence. The maximum deviation among the moments 
obtained with different deuteron wave functions will be taken as an estimate of the 
systematic error due to the wave function model dependence. 

An independent test of the nucleon distribution functions may be made by consider- 
ing quasi-elastic scattering. Within the impulse approximation and neglecting fmite-Q 2 
effects, the quasi-elastic contribution to the deuteron structure function, F% , can 
be written in factorized form, in terms of the nucleon distribution function f D and the 
nucleon elastic form factors, 

' ' ' <15) 

where xd = Q 2 /2Mdv ~ x/2. Thus, one can define phenomenological nucleon distribu- 
tion moments J-^iQ 2 ) in terms of the moments (Q 2 ) of the deuteron quasi-elastic 
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structure function fP^ E ' and the nucleon form factors, 



(l + r)M n D(QE) (Q 2 ) 



G p E \Q 2 ) + G% 2 (Q 2 ) + ri&fm + G^m) 



P 2, 



(16) 



We expect that T^{Q 2 ) ^ T% at large enough Q 2 . Note that for the moment M% (QE) (Q 2 ) 
we use the Nachtmann definition (2H1 to eliminate the finite-Q 2 effects due to the nonzero 
target mass (i.e. the deuteron mass), 



M^ E \Q 2 )^ / dx D ^-F 2 D ( QE \ XD ,Q 2 ) 
Jo x A D 



2 3 + 3(n + l)r D + n(n + 2)r 2 D 



(n + 2)(n + 3) 



(17) 



where r D = Jl + AM 2 D x 2 D /Q 2 and £ D = 2sd/(1 + r D ) 
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Figure 3: Moments (Q 2 ) (a) andF®(Q 2 ) (b) extracted from the quasi- elastic peak moments 
defined in Eq. ^16j) (triangles). The hatched areas represent the corresponding moments and 
T§ calculated from Eq. \1J\) , including the nuclear model uncertainties shown in Fig. 



The deuteron quasi-elastic contribution to the structure function F 2 has been 
calculated using the approach developed in Ref. [21], which includes final state interac- 
tions and has been tested against the recent deuteron experimental data of Ref. As 
an illustration, the extracted J-^iQ 2 ) and J-^(Q 2 ) moments are shown in Fig. |3] as a 
function of Q 2 , and compared with the corresponding moments T 2 and J-§ calculated 
from the distribution function f D (z) in Eq. (|14l) . The results agree well within the 
nuclear model uncertainties. 
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2.3 Off-shell effects 



The expressions for the deuteron structure function thus far have been in the framework 
of the one-dimensional convolution formula, in which nucleon off-shell, relativistic, and 
other effects beyond the impulse approximation are neglected. While this constitutes a 
reasonable approximation at moderate x, or for low moments, at large (and very small) 
x values these effects will become more important. 

Within the nonrelativistic approach discussed in Sec. l2.1,Tl the fact that the structure 
of bound nucleons can differ from that of free nucleons is incorporated in the dependence 
of the bound nucleon structure function in Eq. @ on the nucleon virtuality, or off- 
shellness, p 2 . Since the deuteron is treated as a nonrelativistic system, the off-shellness 
parameter v = (p 2 — M 2 )/M 2 is on average a small number. Expanding the leading 
twist structure function in v and keeping only the first order correction in v, one finds 

F 2 N (x,Q 2 ,p 2 ) = F 2 N (x,Q 2 )(l + 5f 2 (x)v) , (18) 
5f 2 = dlnF^/dlnp 2 , (19) 

where the first term on the right in Eq. (|18|) is the structure function of the on-mass-shell 
nucleon, and the derivative is evaluated at p 2 = M 2 . The deuteron structure function 
F 2 D can then be written as in Eq. ()10j) . with 

8^F 2 D (x,Q 2 ) = + I*d(p)| 2 5f 2 (x/z,Q 2 )vF 2 N (x/z,Q 2 ). (20) 

The off-shell correction Sf 2 was studied phenomenologically in Ref. [23] by analyzing 
data on the ratios of structure functions of different nuclei (EMC ratios). It was found 
that to a good approximation Sf 2 was independent of Q 2 , and could be parameterized 
as 

5f 2 = C N (x - x 1 )(x - x )(l + x - x) , (21) 

with the best fit parameters Cn = 8.1 ± 0.3 s t a t ± 0.5 sys t, x\ = 0.05, and xo = 0.448 ± 
0.005 s tat =t 0.007 sys t. Futhermore, the off-shell correction to the structure function F% is 
given by the same function 5f 2 , which is consistent with the normalization of nuclear 
valence quark distribution (for more details see Ref. 25 ). 

Expressing the off-shell correction to F 2 in terms of moments, one can change the 
order of the p and x integrations to obtain 

M n D (Q 2 ) = jrD M ^(Q 2 ) + 5C n 5M^(Q 2 ) , (22) 

where the first term on the right hand side is the usual convolution term and the second 
is the off-shell correction, with 

ac, » = /(w( 1 + s) 2 "" l,| * D(p)|2 (23) 

and 

5M^(Q 2 ) = f dxx n - 2 Sf 2 (x) F 2 N (x,Q 2 ) . (24) 
J o 
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Writing the nucleon moment in terms of the proton and neutron moments explicitly, 
M£ = (Mn + M™)/2, one can solve Eq. 1)221) to express the neutron moment in terms 
of the experimental proton and deuteron moments, and the theoretical and off-shell 
corrections, 

ilC(Q 2 ) = 2M®(Q 2 ) J " M£(Q 2 ) , (25) 



where the off-shell ratio A^ off ^ is given by 

A "° ff)(g2) = J^I dX ^ 5{0S){X) F ^M 2 ) > (26) 
and 5( off ) is defined as 



From Eqs. (|20j) . (|23j) and (|24|) one can also write the moment correction as 

(28) 



For the relativistic model of Sec, 12.1,21 the off-shell correction S^'F® receives two 
contributions: one which depends on the relativistic P-state wave functions ip± s and ipu, 
and one which is given by the off-shell part of the truncated nucleon structure functions 
Wo, 1,2- The former was found to be of the order 1% and negative over the whole x range, 
while the latter was somewhat smaller in magnitude. The total off-shell ratio 5^°^ was 
parameterized at Q 2 f» 5 GeV 2 as ^3] 

8^(x,Q 2 ) = a (l + a lX a *)(l-(a 3 -x a *D , (29) 

with the parameters ao = —0.014, a\ = 3.0, 02 = 20.0, 03 = 1.067, 04 = 1.5, and 
a?, = 18.0. The Q 2 dependence of this ratio was found to be weak. Note that the 
estimate of 5^ in Ref. [H| was based on the set of wave functions in Ref. derived 
from an NN potential with a pseudoscalar irNN coupling, which are known to give 
larger P-state contributions than more modern treatments |10|ll2j. The correction 5^ oS ^ 
in Eq. Q29|) should therefore serve as an upper limit on the size of the off-shell correction 
in this model. 

The x dependence of the corrections 5^°^ in the models of Refs. |14| I25j . evaluated 
using the leading twist parameterization [Sj of P 2 D a t Q 2 = 10 GeV 2 , is displayed in Fig.0J 
Both of the corrections are around 1-2% for x < 0.6 — 0.7, but increase significantly in 
magnitude for x > 0.8. The correction at large x is negative in both cases, and is 
especially dramatic for the model of Ref. j2E] ■ The resulting moments A^ off ' are listed in 
Tablenf° r re = 2 — 10. The magnitude of the correction for the lowest moment is similar 
in the two models, but is of opposite sign. For higher n, the relative importance of the 
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Figure 4: Nucleon off-shell correction 5^ oS \x,Q 2 ) as a function of x at Q 2 = 10 GeV 2 , in the 
nonrelativistic fiffl (dashed) and relativistic \14\ (solid) models. 



Table 1: Moments of the non- convolution corrections to the deuteron structure function 
for n = 2 — 10: the first two columns give off-shell corrections Al° S ^ in the relativistic \14\ and 
nonrelativistic ]25$ models, respectively; the third and fourth columns represent the contributions 
from meson exchange currents and nuclear shadowing respectively. 



n 


A { n m 14 


A^f } 25, 




A (shad) 


2 


-5.67xl0" 3 


5.08X10" 3 


2.6 xitr 4 


-1.0 xlO" 4 


4 


-1.21xl0~ 2 


-3.34xl0~ 3 


7.6 xio~ e 


7.6 xl0~ 6 


6 


-1.39xHT 2 


-1.68xl0~ 2 


8.4 xKT 7 


1.5 xl0~ s 


8 


-1.46x10^ 


-3.01xl0" 2 


1.6 xlO" 7 


4.1 xlO" 7 


10 


-1.53xl0 -2 


-3.97xl0~ 2 


4.4 xl(T s 


1.4xl0 -7 



off-shell correction increases. Because of the steep decrease of 5^ oS ^ in the nonrelativistic 
model |25j . the magnitude of the correction An^ at n = 10 is about 2-3 times larger 
than the one in the model of Ref . |14j . 

Note that part of the off-shell corrections in the model of Ref. ^1] are included in 
the relativistic distribution function (|llj). which in the model of Ref. are contained 
only in the non-convolution term. Therefore the definition of off-shell corrections can 
differ between various models, and one should ensure that off-shell effects in general are 
calculated self-consistently. The difference between the two models provides an estimate 
of the systematic error associated with the poor knowledge of the change in the nucleon 
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structure when the nucleon is bound in deuteron (see also Ref. [26 ). 

2.4 Small- a; effects: shadowing and meson exchange currents 

For completeness, in this section we briefly discuss the nuclear corrections in the deuteron 
at small values of x. From the definition of the moments of the F2 structure function, it 
is clear that small- x effects will mostly be relevant for the lowest moment, and will be 
negligible for n > 4. 

At small values of x, the coherence length of a virtual qq fluctuation of the photon 
probe in the laboratory frame is ~ 1/Mx, which for x < 0.1 becomes larger than the 
inter-nucleon separation. Consequently at small x the scattering from the deuteron can 
take place via double scattering from both nucleons, which gives rise to the phenomenon 
of nuclear shadowing (see e.g. Refs. [UE1)- 

At low Q 2 the qq fluctuations are usually represented in the form of the vector mesons, 
p, to and <f>. Their contributions scale as 1/Q 2 , however, and are formally of higher twist. 
Leading twist shadowing contributions arise from diffractive scattering of the virtual 
photon from the bound nucleons, which is described through Pomeron exchange, as well 
as contributions in which the exchanged quanta are mesons (most notably the pion). 

0.002 
0.001 





u -0.001 

1 

To -0.002 
-0.003 
-0.004 



0.1 0.2 0.3 0.4 0.5 

x 

Figure 5: Shadowing (Pomeron exchange) and meson exchange current corrections S^ p \x,Q 2 ) 
(solid) and <$( MEC ) (x, Q 2 ) (dotted) as a function of x at Q 2 = 10 GeV 2 0/. 

Both the Pomeron (P) exchange and meson exchange current (MEC) contributions 
to the deuteron structure function have been estimated by several authors P3 El El 
0JEH1EI!) and enter as additive corrections to FrP, 

F 2 D (x,Q 2 ) - F? (z,Q 2 ) + S( p ^F 2 D (*,Q 2 ) . (30) 
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The Pomeron shadowing contribution S^F® is negative, while the meson (pion) ex- 
change contribution S^ MEC ^F^ is positive and hence corresponds to an antishadowing 
effect. To estimate their effect on the moment analysis, we use the results from Ref. 
which can be conveniently parameterized at Q 2 ~ 5 GeV 2 as 

j( p )F 2 D = b x bl (l - x) b2 (l + b 3 x b *) , (31) 

for the P-exchange, where bo = —0.003, b\ = —0.13, 62 = 5.0, ^3 = —2.2 and 64 = 0.4, 
and 

6 (MEC) F D = CqX c 1{1 _ x) c 2 (32) 

for meson exchange, with cq = 0.002, c\ = 0.03 and C2 = 6.0. As a fraction of the 
empirical leading twist F®, the ratios 

«W» >(x,Q') = ^|^2 (33) 
are plotted in Fig. [5J The corresponding moments, defined in analogy with Eq. (|26|). 

a (p,mec )(q2) = ^fdx x- 2 ^* MEC )(z, Q 2 ) F 2 D (x, Q 2 ) , (34) 
are listed in Tabled where it is clear that these corrections are negligible for n > 4. 

3 Moments of the neutron structure function 

In this section we present the extraction of the lowest few moments of the neutron 
structure function F£, taking into account the systematic uncertainties associated with 
our incomplete knowledge of nuclear effects in the deuteron. We shall focus on the leading 
twist (LT) part of the structure function, and make use of the previous determination 
of the LT moments of the proton and deuteron structure functions obtained by the 
twist analyses of Refs. Higher twist contributions to the moments, in the form 

of subleading 1/Q 2 power corrections, will be discussed in Sec. El Moreover, since the 
twist extraction performed in Refs. |5J |Sj is based on the Nachtmann definition of the 
moments j2Hj, our LT terms are determined without any a priori knowledge of the re- 
shape of parton distribution functions (PDFs). 

The moments of the LT F% structure function of the neutron are related to those of 
proton and deuteron by 

M-(Q 2 ) = 2M°(Q 2 ) 1 {Q) TD An - M£(Q 2 ) , (35) 

where !F® /[l — Ai° ff ^(Q 2 ) — A^' MEC \q 2 )] represents the nuclear correction factor, or in 
other words the EMC effect of the deuteron in moment space. From the values reported 
in Fig. 12 and in Table Q the EMC effect on the second moment is expected to be below 
1% [see also Fig.ES>)]. 
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Using the nuclear ingredients discussed in the previous Section, the extracted LT 
moments M™(Q 2 ) are listed in Table El and shown in Fig. El The central values refer 
to the case of the relativistic distribution Hll[) with the deuteron wave function from 
Ref. |12j . For comparison we also plot the corresponding proton (stars) and deuteron 
(open crosses) LT moments. It turns out that the neutron LT moments with n < 8 
can be determined with quite good precision, namely < 18% (statistical) and < 30% 
(systematic). Different contributions to the total systematic error are shown in Fig. [7| 
It can be seen that all the uncertainties related to the nuclear correction factor steeply 
rise with n, so that the total errors due to nuclear effects become comparable with the 
experimental ones around n = 10. Thus, an improvement in the measurements of the 
proton and deuteron structure functions at large x would be beneficial and also the 
moments at n = 12 may be still reasonably accessible. 

In Fig. |H1 our extracted neutron moments are compared with the corresponding mo- 
ments calculated from various PDF sets available in the literature |281 1291 1301 13 lj . For 
n > 2 all PDF sets predict significantly lower moments at low Q 2 . This difference is 
due to Sudakov effects at large-x; indeed, as was firstly investigated in Ref. 32_, soft 
gluon resummation plays a relevant role for n > 2 already at Q 2 = 3 — 4 GeV 2 . While 
resummation effects are taken into account in the analysis of Refs. 0101) a h the PDF 
sets of Refs. [2H1 EH1 E01 EU are based on next-to-leading (NLO) or NNLO treatments. 
At higher Q 2 the PDF predictions are consistent with our results for n < 10 within the 
quoted statistical errors. For M2 however all PDFs give somehow smaller values at large 
Q 2 , though within the statistical uncertainties. 

The large-n behavior of the moments is also closely related to the asymptotic x — > 1 
behavior of the structure functions. Indeed, inspecting the definition of the moments, 
Eq. @, we conclude that at n > 1 the factor x n ~ 2 in the integrand suppresses the 
contribution from the region of small x, leaving the moment effectively saturated by a 
region close to x = 1. Furthermore, assuming that the structure functions F% and F^ 
vanish as (1 — x)@ with the same exponent for the proton and neutron (as one would 
expect from perturbative QCD) [331, it is possible to derive the asymptotic relation 
between the ratios of the moments and the structure functions 

lim — % = lim . 36 

This relation can be useful in obtaining some insights into the large- x asymptotics of the 
d/u ratio from the large-n limit of the neutron and proton moments. Note that most of 
the PDF fits [SHlEni EDI EU are characterized by the limiting value lim^-i F$/F$ = 1/4, 
which is expected if soft, nonperturbative effects dominate the scattering at large x 
(namely, if the energy of the spectator uu system which accompanies the d quark is 
increased relative to the ud system which is relevant for the u quark [25)- O n the 
other hand, arguments based on helicity conservation in pQCD, when combined with an 
unperturbed spin-flavor symmetric wave function, suggest an asymptotic value of 3/7 

m 

The ratios of neutron to proton LT moments are shown in Fig.|21for several n, where 
they are also compared to various PDF parameterizations |28l 1291 1301 131j . The ratio 
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Table 2: Extracted leading twist moments of the neutron F% structure functions, M™(Q 2 ), for 
n = 2,4,6,8, shown in Fig. together with statistical and systematic uncertainties. 



— TV? — r7n — 771 — 771 

Q [(GeV/c)^J 


M£(Q' Z ) x 10 _1 




M£(Q' Z ) x 10~ d 


M£(Q' Z ) x 10~ d 


1.025 


1.37 ± 0.06 ± 0.12 


1.95 zh 0.20 ± 0.11 


8.06 ± 0.45 ± 1.15 


5.51 zh 0.98 zh 1.49 


1.075 


1.36 ± 0.06 ± 0.12 


1.88 ± 0.19 ± 0.11 


7.35 ± 0.41 ± 1.04 


4.59 zh 0.81 zh 1.24 


1.125 


1.36 ± 0.06 ± 0.12 


1.82 ± 0.18 ± 0.11 


6.78 ± 0.38 ± 0.96 


3.95 zh 0.70 zh 1.07 


1.175 


1.36 zh 0.06 ± 0.12 


1.76 zh 0.18 ± 0.10 


6.33 zh 0.35 zh 0.90 


3.48 zh 0.62 zh 0.94 


1.225 


1.36 zh 0.05 ± 0.11 


1.71 zh 0.17 ± 0.10 


5.95 ± 0.33 ± 0.85 


3.13 zh 0.56 zh 0.85 


1.275 


1.35 zh 0.05 ± 0.11 


1.67 zh 0.17 ± 0.10 


5.63 ± 0.31 ± 0.80 


2.85 zh 0.51 zh 0.77 


1.325 


1.35 ± 0.05 ± 0.11 


1.63 ± 0.16 ± 0.09 


5.36 zh 0.30 zh 0.76 


2.63 zh 0.47 zh 0.71 


1.375 


1.35 ± 0.05 ± 0.11 


1.59 ± 0.16 ± 0.09 


5.13 ± 0.28 zh 0.73 


2.44 zh 0.43 zh 0.66 


1.425 


1.34 ± 0.05 ± 0.11 


1.56 ± 0.16 ± 0.09 


4.93 ± 0.27 ± 0.70 


2.29 zh 0.41 zh 0.62 


1.475 


1.34 zh 0.05 ± 0.11 


1.53 ± 0.15 ± 0.09 


4.74 zh 0.26 zh 0.68 


2.16 zh 0.38 zh 0.58 


1.525 


1.34 zh 0.05 zh 0.11 


1.50 zh 0.15 zh 0.09 


4.58 zh 0.25 zh 0.65 


2.04 zh 0.36 zh 0.55 


1.575 


1.34 zh 0.05 ± 0.11 


1.48 zh 0.15 ± 0.09 


4.44 zh 0.25 zh 0.63 


1.94 zh 0.35 zh 0.53 


1.625 


1.33 zh 0.05 ± 0.11 


1.45 ± 0.15 ± 0.08 


4.31 ± 0.24 ± 0.61 


1.86 zh 0.33 zh 0.50 


1.675 


1.33 ± 0.05 ± 0.11 


1.43 ± 0.14 ± 0.08 


4.19 ± 0.23 ± 0.60 


1.78 zh 0.32 zh 0.48 


1.725 


1.33 ± 0.05 ± 0.11 


1.41 ± 0.14 ± 0.08 


4.08 zh 0.23 zh 0.58 


1.71 zh 0.30 zh 0.46 


1.775 


1.32 ± 0.05 ± 0.11 


1.39 ± 0.14 ± 0.08 


3.99 ± 0.22 ± 0.57 


1.65 zh 0.29 zh 0.45 


1.825 


1.32 ± 0.05 ± 0.11 


1.37 ± 0.14 ± 0.08 


3.90 ± 0.22 ± 0.55 


1.60 zh 0.28 zh 0.43 


1.875 


1.32 zh 0.05 ± 0.11 


1.36 zh 0.14 ± 0.08 


3.81 zh 0.21 ± 0.54 


1.55 zh 0.27 zh 0.42 


1.925 


1.32 zh 0.05 ± 0.11 


1.34 ± 0.13 ± 0.08 


3.73 ± 0.21 ± 0.53 


1.50 zh 0.27 zh 0.41 


1.975 


1.32 ± 0.05 ± 0.11 


1.33 ± 0.13 ± 0.08 


3.66 ± 0.20 ± 0.52 


1.46 zh 0.26 zh 0.40 


2.025 


1.31 zh 0.05 ± 0.11 


1.31 zh 0.13 ± 0.08 


3.59 zh 0.20 zh 0.51 


1.42 zh 0.25 zh 0.38 


2.075 


1.31 zh 0.05 ± 0.11 


1.30 ± 0.13 ± 0.08 


3.53 ± 0.20 ± 0.50 


1.38 zh 0.25 zh 0.38 


2.125 


1.31 ± 0.05 ± 0.11 


1.29 ± 0.13 ± 0.07 


3.47 ± 0.19 ± 0.49 


1.35 zh 0.24 zh 0.37 


2.175 


1.31 zh 0.05 ± 0.10 


1.27 ± 0.13 ± 0.07 


3.42 ± 0.19 ± 0.49 


1.32 zh 0.23 zh 0.36 


2.225 


1.31 zh 0.05 ± 0.10 


1.26 zh 0.13 ± 0.07 


3.36 zh 0.19 ± 0.48 


1.29 zh 0.23 zh 0.35 


2.275 


1.31 zh 0.05 ± 0.10 


1.25 zh 0.13 ± 0.07 


3.32 zh 0.18 ± 0.47 


1.27 zh 0.22 zh 0.34 


2.325 


1.31 ± 0.05 ± 0.10 


1.24 ± 0.12 ± 0.07 


3.28 zh 0.18 ± 0.47 


1.25 zh 0.22 zh 0.34 


2.375 


1.30 zh 0.05 ± 0.10 


1.23 ± 0.12 ± 0.07 


3.24 ± 0.18 ± 0.46 


1.23 zh 0.22 zh 0.33 


2.425 


1.30 zh 0.05 ± 0.10 


1.22 ± 0.12 ± 0.07 


3.20 ± 0.18 ± 0.46 


1.21 zh 0.21 zh 0.33 


2.475 


1.30 ± 0.05 ± 0.10 


1.22 ± 0.12 ± 0.07 


3.17 ± 0.18 ± 0.45 


1.19 zh 0.21 zh 0.32 


2.525 


1.30 ± 0.05 ± 0.10 


1.21 ± 0.12 ± 0.07 


3.13 ± 0.17 zh 0.45 


1.17 zh 0.21 zh 0.32 


2.575 


1.30 zh 0.05 ± 0.10 


1.20 ± 0.12 ± 0.07 


3.10 zh 0.17 ± 0.44 


1.15 zh 0.21 zh 0.31 


2.625 


1.30 zh 0.05 ± 0.10 


1.19 zh 0.12 zh 0.07 


3.07 zh 0.17 dz 0.44 


1.14 zh 0.20 zh 0.31 


2.675 


1.30 zh 0.05 ± 0.10 


1.18 zh 0.12 ± 0.07 


3.04 ± 0.17 zh 0.43 


1.12 zh 0.20 zh 0.30 


2.725 


1.30 zh 0.05 ± 0.10 


1.18 ± 0.12 ± 0.07 


3.01 zh 0.17 zh 0.43 


1.11 zh 0.20 zh 0.30 


2.775 


1.30 zh 0.05 ± 0.10 


1.17 zh 0.12 ± 0.07 


2.98 zh 0.17 zh 0.42 


1.10 zh 0.19 zh 0.30 


2.825 


1.30 ± 0.05 ± 0.10 


1.16 ± 0.12 ± 0.07 


2.96 zh 0.16 zh 0.42 


1.08 zh 0.19 zh 0.29 


2.875 


1.30 ± 0.05 ± 0.10 


1.16 ± 0.12 ± 0.07 


2.93 zh 0.16 zh 0.42 


1.07 zh 0.19 zh 0.29 


2.925 


1.30 zh 0.05 ± 0.10 


1.15 ± 0.12 ± 0.07 


2.91 zh 0.16 zh 0.41 


1.06 zh 0.19 zh 0.29 


2.975 


1.30 zh 0.05 ± 0.10 


1.15 zh 0.12 ± 0.07 


2.88 zh 0.16 zh 0.41 


1.05 zh 0.19 zh 0.28 


3.025 


1.30 zh 0.05 ± 0.10 


1.14 ± 0.11 ± 0.07 


2.86 zh 0.16 zh 0.41 


1.04 zh 0.18 zh 0.28 


3.075 


1.30 zh 0.05 ± 0.10 


1.13 zh 0.11 ± 0.07 


2.84 zh 0.16 zh 0.40 


1.03 zh 0.18 zh 0.28 


3.125 


1.30 zh 0.05 ± 0.10 


1.13 ± 0.11 ± 0.07 


2.82 zh 0.16 zh 0.40 


1.02 zh 0.18 zh 0.28 


3.175 


1.30 ± 0.05 ± 0.10 


1.12 ± 0.11 ± 0.07 


2.80 zh 0.16 zh 0.40 


1.01 zh 0.18 zh 0.27 


3.225 


1.30 ± 0.05 ± 0.10 


1.12 ± 0.11 ± 0.06 


2.78 zh 0.15 zh 0.40 


1.00 zh 0.18 zh 0.27 


3.275 


1.30 zh 0.05 ± 0.10 


1.11 zh 0.11 ± 0.06 


2.76 zh 0.15 zh 0.39 


0.99 zh 0.18 zh 0.27 


3.325 


1.30 zh 0.05 ± 0.10 


1.11 zh 0.11 ± 0.06 


2.74 zh 0.15 zh 0.39 


0.98 zh 0.17 zh 0.27 


3.375 


1.30 ± 0.05 ± 0.10 


1.10 ± 0.11 ± 0.06 


2.72 ± 0.15 ± 0.39 


0.97 ± 0.17 ± 0.26 


3.425 


1.30 zh 0.05 ± 0.10 


1.10 ± 0.11 ± 0.06 


2.70 zh 0.15 zh 0.38 


0.96 zh 0.17 zh 0.26 


3.475 


1.30 ± 0.05 ± 0.10 


1.09 ± 0.11 ± 0.06 


2.68 zh 0.15 zh 0.38 


0.95 zh 0.17 zh 0.26 


3.525 


1.29 ± 0.05 ± 0.10 


1.09 ± 0.11 ± 0.06 


2.67 zh 0.15 zh 0.38 


0.94 zh 0.17 zh 0.26 


3.575 


1.29 ± 0.05 ± 0.10 


1.08 ± 0.11 ± 0.06 


2.65 zh 0.15 zh 0.38 


0.94 zh 0.17 zh 0.25 


3.625 


1.29 ± 0.05 ± 0.10 


1.08 zh 0.11 ± 0.06 


2.64 zh 0.15 zh 0.38 


0.93 zh 0.17 zh 0.25 


3.675 


1.29 ± 0.05 ± 0.10 


1.08 zh 0.11 ± 0.06 


2.62 zh 0.15 zh 0.37 


0.92 zh 0.16 zh 0.25 


3.725 


1.29 zh 0.05 ± 0.10 


1.07 ± 0.11 ± 0.06 


2.60 zh 0.14 zh 0.37 


0.92 zh 0.16 zh 0.25 


3.775 


1.29 ± 0.05 ± 0.10 


1.07 zh 0.11 ± 0.06 


2.59 zh 0.14 zh 0.37 


0.91 zh 0.16 zh 0.25 


3.825 


1.29 ± 0.05 ± 0.10 


1.06 ± 0.11 ± 0.06 


2.58 zh 0.14 zh 0.37 


0.90 zh 0.16 zh 0.24 
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Q 2 GcV/c y] 


M%(Q' Z ) x 10 _1 


M"(Q^) x 10 - ^ 


MV-(Q' Z ) x lO - ' 3 


Mg{Q' z ) x 10~- 1 


3.875 


1.29 ± 0.05 ± 0.10 


1.06 ± 0.11 ± 0.06 


2.56 ± 0.14 ± 0.36 


0.90 ± 0.16 ± 0.24 


3.925 


1.29 ± 0.05 ± 0.10 


1.06 ± 0.11 ± 0.06 


2.55 ± 0.14 ± 0.36 


0.89 ± 0.16 ± 0.24 


3.975 


1.29 ± 0.05 ± 0.10 


1.05 ± 0.11 ± 0.06 


2.53 ± 0.14 ± 0.36 


0.88 ± 0.16 ± 0.24 


4.025 


1.29 ± 0.05 ± 0.10 


1.05 ± 0.11 ± 0.06 


2.52 ± 0.14 ± 0.36 


0.88 ± 0.16 ± 0.24 


4.075 


1.29 ± 0.05 ± 0.10 


1.05 ± 0.11 ± 0.06 


2.51 ± 0.14 ± 0.36 


0.87 ± 0.15 ± 0.24 


4.125 


1.29 ± 0.05 ± 0.10 


1.04 ± 0.10 ± 0.06 


2.50 ± 0.14 ± 0.36 


0.87 ± 0.15 ± 0.23 


4.175 


1.29 ± 0.05 ± 0.10 


1.04 ± 0.10 ± 0.06 


2.48 ± 0.14 ± 0.35 


0.86 ± 0.15 ± 0.23 


4.225 


1.29 ± 0.05 ± 0.10 


1.04 ± 0.10 ± 0.06 


2.47 ± 0.14 ± 0.35 


0.86 ± 0.15 ± 0.23 


4.275 


1.29 ± 0.05 ± 0.10 


1.03 ± 0.10 ± 0.06 


2.46 ± 0.14 ± 0.35 


0.85 ± 0.15 ± 0.23 


4.325 


1.29 ± 0.05 ± 0.10 


1.03 ± 0.10 ± 0.06 


2.45 ± 0.14 ± 0.35 


0.84 ± 0.15 ± 0.23 


4.375 


1.29 ± 0.05 ± 0.10 


1.03 ± 0.10 ± 0.06 


2.44 ± 0.14 ± 0.35 


0.84 ± 0.15 ± 0.23 


4.425 


1.29 ± 0.05 ± 0.10 


1.02 ± 0.10 ± 0.06 


2.43 ± 0.13 ± 0.35 


0.84 ± 0.15 ± 0.23 


4.475 


1.29 ± 0.05 ± 0.10 


1.02 ± 0.10 ± 0.06 


2.42 ± 0.13 ± 0.34 


0.83 ± 0.15 ± 0.23 


4.525 


1.29 ± 0.05 ± 0.10 


1.02 ± 0.10 ± 0.06 


2.41 ± 0.13 ± 0.34 


0.83 ± 0.15 ± 0.22 


4.575 


1.29 ± 0.05 ± 0.10 


1.01 ± 0.10 ± 0.06 


2.40 ± 0.13 ± 0.34 


0.82 ± 0.15 ± 0.22 


4.625 


1.29 ± 0.05 ± 0.10 


1.01 ± 0.10 ± 0.06 


2.39 ± 0.13 ± 0.34 


0.82 ± 0.15 ± 0.22 


4.675 


1.29 ± 0.05 ± 0.10 


1.01 ± 0.10 ± 0.06 


2.38 ± 0.13 ± 0.34 


0.81 ± 0.14 ± 0.22 


4.725 


1.29 ± 0.05 ± 0.10 


1.01 ± 0.10 ± 0.06 


2.37 ± 0.13 ± 0.34 


0.81 ± 0.14 ± 0.22 


4.775 


1.29 ± 0.05 ± 0.10 


1.00 ± 0.10 ± 0.06 


2.36 ± 0.13 ± 0.34 


0.80 ± 0.14 ± 0.22 


4.825 


1.29 ± 0.05 ± 0.10 


1.00 ± 0.10 ± 0.06 


2.35 ± 0.13 ± 0.33 


0.80 ± 0.14 ± 0.22 


4.875 


1.29 ± 0.05 ± 0.10 


1.00 ± 0.10 ± 0.06 


2.34 ± 0.13 ± 0.33 


0.80 ± 0.14 ± 0.22 


4.925 


1.29 ± 0.05 ± 0.10 


1.00 ± 0.10 ± 0.06 


2.33 ± 0.13 ± 0.33 


0.79 ± 0.14 ± 0.21 


4.975 


1.29 ± 0.05 ± 0.10 


0.99 ± 0.10 ± 0.06 


2.32 ± 0.13 ± 0.33 


0.79 ± 0.14 ± 0.21 


5.025 


1.29 ± 0.05 ± 0.10 


0.99 ± 0.10 ± 0.06 


2.31 ± 0.13 ± 0.33 


0.78 ± 0.14 ± 0.21 


5.075 


1.29 ± 0.05 ± 0.10 


0.99 ± 0.10 ± 0.06 


2.30 ± 0.13 ± 0.33 


0.78 ± 0.14 ± 0.21 


5.125 


1.29 ± 0.05 ± 0.10 


0.99 ± 0.10 ± 0.06 


2.30 ± 0.13 ± 0.33 


0.78 ± 0.14 ± 0.21 


5.275 


1.29 ± 0.05 ± 0.10 


0.98 ± 0.10 ± 0.06 


2.27 ± 0.13 ± 0.32 


0.77 ± 0.14 ± 0.21 


5.325 


1.29 ± 0.05 ± 0.10 


0.98 ± 0.10 ± 0.06 


2.26 ± 0.13 ± 0.32 


0.76 ± 0.14 ± 0.21 


5.375 


1.29 ± 0.05 ± 0.10 


0.97 ± 0.10 ± 0.06 


2.26 ± 0.13 ± 0.32 


0.76 ± 0.13 ± 0.21 


5.475 


1.28 ± 0.05 ± 0.09 


0.97 ± 0.10 ± 0.06 


2.24 ± 0.12 ± 0.32 


0.75 ± 0.13 ± 0.20 


5.525 


1.28 ± 0.05 ± 0.09 


0.97 ± 0.10 ± 0.06 


2.23 ± 0.12 ± 0.32 


0.75 ± 0.13 ± 0.20 


5.625 


1.28 ± 0.05 ± 0.09 


0.96 ± 0.10 ± 0.06 


2.22 ± 0.12 ± 0.32 


0.74 ± 0.13 ± 0.20 


5.675 


1.28 ± 0.05 ± 0.09 


0.96 ± 0.10 ± 0.06 


2.21 ± 0.12 ± 0.31 


0.74 ± 0.13 ± 0.20 


5.725 


1.28 ± 0.05 ± 0.09 


0.96 ± 0.10 ± 0.06 


2.21 ± 0.12 ± 0.31 


0.74 ± 0.13 ± 0.20 


5.955 


1.28 ± 0.05 ± 0.09 


0.95 ± 0.10 ± 0.06 


2.18 ± 0.12 ± 0.31 


0.72 ± 0.13 ± 0.20 


6.915 


1.28 ± 0.04 ± 0.09 


0.92 ± 0.09 ± 0.05 


2.07 ± 0.12 ± 0.29 


0.68 ± 0.12 ± 0.18 


7.267 


1.28 ± 0.04 ± 0.09 


0.91 ± 0.09 ± 0.05 


2.04 ± 0.11 ± 0.29 


0.67 ± 0.12 ± 0.18 


7.630 


1.28 ± 0.04 ± 0.09 


0.90 ± 0.09 ± 0.05 


2.00 ± 0.11 ± 0.29 


0.65 ± 0.12 ± 0.18 


8.021 


1.28 ± 0.04 ± 0.09 


0.89 ± 0.09 ± 0.05 


1.97 ± 0.11 ± 0.28 


0.64 ± 0.11 ± 0.17 


8.847 


1.28 ± 0.04 ± 0.09 


0.87 ± 0.09 ± 0.05 


1.91 ± 0.11 ± 0.27 


0.62 ± 0.11 ± 0.17 


9.775 


1.27 ± 0.04 ± 0.09 


0.85 ± 0.09 ± 0.05 


1.86 ± 0.10 ± 0.26 


0.59 ± 0.11 ± 0.16 


10.267 


1.27 ± 0.04 ± 0.09 


0.85 ± 0.09 ± 0.05 


1.83 ± 0.10 ± 0.26 


0.58 ± 0.10 ± 0.16 


10.762 


1.27 ± 0.04 ± 0.09 


0.84 ± 0.08 ± 0.05 


1.81 ± 0.10 ± 0.26 


0.57 ± 0.10 ± 0.16 


11.344 


1.27 ± 0.04 ± 0.09 


0.83 ± 0.08 ± 0.05 


1.78 ± 0.10 ± 0.25 


0.56 ± 0.10 ± 0.15 


12.580 


1.27 ± 0.04 ± 0.09 


0.81 ± 0.08 ± 0.05 


1.73 ± 0.10 ± 0.25 


0.54 ± 0.10 ± 0.15 


13.238 


1.27 ± 0.04 ± 0.09 


0.80 ± 0.08 ± 0.05 


1.71 ± 0.10 ± 0.24 


0.53 ± 0.09 ± 0.15 


14.689 


1.27 ± 0.04 ± 0.09 


0.79 ± 0.08 ± 0.05 


1.66 ± 0.09 ± 0.24 


0.52 ± 0.09 ± 0.14 


17.108 


1.26 ± 0.04 ± 0.09 


0.77 ± 0.08 ± 0.04 


1.60 ± 0.09 ± 0.23 


0.49 ± 0.09 ± 0.13 


19.072 


1.26 ± 0.04 ± 0.09 


0.75 ± 0.08 ± 0.04 


1.56 ± 0.09 ± 0.22 


0.48 ± 0.08 ± 0.13 


20.108 


1.26 ± 0.04 ± 0.09 


0.75 ± 0.08 ± 0.04 


1.54 ± 0.09 ± 0.22 


0.47 ± 0.08 ± 0.13 


21.097 


1.26 ± 0.04 ± 0.09 


0.74 ± 0.07 ± 0.04 


1.52 ± 0.08 ± 0.22 


0.46 ± 0.08 ± 0.13 


24.259 


1.26 ± 0.04 ± 0.08 


0.73 ± 0.07 ± 0.04 


1.47 ± 0.08 ± 0.21 


0.44 ± 0.08 ± 0.12 


26.680 


1.26 ± 0.04 ± 0.08 


0.71 ± 0.07 ± 0.04 


1.44 ± 0.08 ± 0.21 


0.43 ± 0.08 ± 0.12 


32.500 


1.25 ± 0.04 ± 0.08 


0.69 ± 0.07 ± 0.04 


1.39 ± 0.08 ± 0.20 


0.41 ± 0.07 ± 0.11 


34.932 


1.25 ± 0.04 ± 0.08 


0.69 ± 0.07 ± 0.04 


1.37 ± 0.08 ± 0.19 


0.41 ± 0.07 ± 0.11 


36.750 


1.25 ± 0.04 ± 0.08 


0.68 ± 0.07 ± 0.04 


1.35 ± 0.08 ± 0.19 


0.40 ± 0.07 ± 0.11 


43.970 


1.25 ± 0.04 ± 0.08 


0.66 ± 0.07 ± 0.04 


1.31 ± 0.07 ± 0.19 


0.38 ± 0.07 ± 0.10 


47.440 


1.25 ± 0.04 ± 0.08 


0.66 ± 0.07 ± 0.04 


1.29 ± 0.07 ± 0.18 


0.38 ± 0.07 ± 0.10 


64.270 


1.25 ± 0.04 ± 0.08 


0.63 ± 0.06 ± 0.04 


1.22 ± 0.07 ± 0.17 


0.35 ± 0.06 ± 0.10 


75.000 


1.25 ± 0.04 ± 0.08 


0.62 ± 0.06 ± 0.04 


1.19 ± 0.07 ± 0.17 


0.34 ± 0.06 ± 0.09 


86.000 


1.24 ± 0.04 ± 0.08 


0.61 ± 0.06 ± 0.04 


1.16 ± 0.06 ± 0.17 


0.33 ± 0.06 ± 0.09 


97.690 


1.24 ± 0.04 ± 0.08 


0.60 ± 0.06 ± 0.03 


1.13 ± 0.06 ± 0.16 


0.32 ± 0.06 ± 0.09 
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Figure 6: Leading twist moments of the neutron structure function for n = 2,4,6 and 8 as 
a function of Q 2 (triangles), extracted using the relativistic distribution in Eq. j 11}) with the 
deuteron wave function from Ref. \lty - For comparison we also show the proton (stars) and 
deuteron (open crosses) moments. The errors shown are statistical. 

depends only slightly on Q 2 for n = 2, while it is almost completely Q 2 -independent for 
larger n. The latter are in reasonable agreement with our extracted ratios within both 
statistical and systematical errors. 

The n/p ratio systematically decreases with n and reaches the value ~ 0.34±0.12 stat ± 
0.13 syst at n — 10. To compare this with the theoretical predictions the mean value of x 
probed by our moments needs to be estimated. This can be done in the following way: 
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Figure 7: Various contributions to the systematic errors on the LT neutron moments for 
n = 2 — 10: deuteron LT moment error (solid); proton LT moment error (dashed); error due to 
the momentum cut (dotted); dependence on the relativistic treatment (dot-dashed); offshellness 
uncertainty (left-hatched); wave function model dependence (cross-hatched) . 



for any calculated n-th neutron and proton moment we estimate a corresponding average 
value of x defined as (x) n = (M™ +1 + M% +1 )/(M% + M%); then the ratio of neutron to 
proton structure functions is computed at the values (x) n . 

The results of this procedure are shown in Fig. (dashed curve), where they are 
also compared with the ratio of the extracted LT moments (solid curve). The two ratios 
coincide with good accuracy, suggesting that the extracted ratio of moments is similar to 
the ratio of structure functions taken at the average value {x) n . Note that at n = 10 one 
has (a) io ^ 0.70, and the F 2 n /Ff ratio at x ~ 0.70 is equal to 0.34 ± 0.12 stat ± 0.13 syst . 
As discussed above (see Fig. |7|, moments for n = 12 may still be accessible by improving 
measurements of the proton and deuteron structure functions at large x, and with n = 12 
one can probe values of x as large as 0.75. However, higher (x) n values are difficult to 
access by this method because of the rapid increase of the uncertainties on the nuclear 
corrections with the order n. 

The sensitivity of the moment's ratio to the high-x tail of the PDFs is illustrated in 
Fig. 1111 where our n/p moment's ratio (triangles) is compared with the PDF parame- 
terizations from Ref. [30] with and without modifying the <i-quark distribution as: 

d(x) ->d(x) + 0.1 x m (l + x) u(x). (37) 

with m a parameter. Depending on the value of m, the form in Eq. (|37|) allows the x — ► 1 
behavior of the d/u ratio to vary between and 1/5, which correspond respectively to 
the structure function limits lirn^i F^ix^jF^ ( x ) °f an d 3/7. The dotted line in 
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0.15 




Q 2 (GeV/c) 2 Q 2 (GeV/c) 2 

Figure 8: Extracted moments of the neutron structure function for n = 2,4,6 and 8 as a 
function of Q 2 in comparison to moments of different PDFs: solid line represents GRV ' L 28l, 
dashed line represents Alekhin dotted line represents CTEQ \S0f . dot-dashed line represents 
MRST \29f . The hatched area shows systematic errors. 



Fig. ^2 corresponds to the case m = 1 in Eq. I|37|). which was advocated in the NNLO 
analysis of both electron and neutrino DIS data of Ref. |36| . It appears that our results 
rule out such an enhancement of the d quark distribution, since the correction is not 
negligible even at values of x as low as x ~ 0.5. A larger value of the parameter m 
confines the enhancement only to larger values of x, and the results obtained adopting 
value m = 4, corresponding to the dashed line in Fig. 111! are more consistent with our 
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Figure 9: Ratio of neutron to proton leading twist moments for various n (triangles), together 
with the statistical and systematic uncertainties (hatched area). The horizontal lines denote the 
ratios obtained with different PDF sets: GRV J^Ej (solid), Alekhin JEjJj (dashed), CTEQ 
(dotted), MRST (dot-dashed). 



extracted n/p ratio. 

Before closing this section, our results for the LT moments of the nonsinglet p — 
n part of the F 2 structure function, Ml~ n {Q 2 ) = M%(Q 2 ) - M"(Q 2 ), are listed in 
Table 01 together with their statistical and systematic uncertainties. We point out that 
an important fraction of the total systematic error for the second moments M^iQ 2 ) and 
M^iQ 2 ) comes from the extrapolation in the unmeasured low-x regions. Such systematic 
uncertainty basically cancels out in the isovector p — n combination. This allows us to 
determine M2~ n (Q 2 ) within ~ 20% global accuracy. 

Low-order moments of the isovector u — d quark distributions have been recently 
calculated in both unquenched [37] and quenched [SHI lattice QCD. Our results for 
the nonsinglet LT moments Mn n can be easily compared with those of the isovector 
moments ((x n ~ 1 ) u — (x n " 1 )^), taking into account that for even values of the order n 
one has 

Mr n (Q 2 ) = \C™ {{x n ~ l ) u - {x n - l ) d ) , (38) 

where C^ s is the (n — l)-th moment of the quark coefficient function. At the scale 
Q 2 = 4 GeV 2 and in the MS renormalization scheme at NLO, which are the scale and 
the scheme adopted in Refs. (33 EB], one gets C^ s ~ 1.01 and Cf s ~ 1.15. 

The quark masses accessible in present lattice calculations are still relatively high and 
thus the extrapolation to the physical point is one of the open issues in lattice applications 
to light hadron phenomenology. In Refs. |37U38| the lattice data are extrapolated linearly 
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Table 3: Leading twist moments of the nonsinglet p — n part of the F 2 structure function, 
MP~ n (Q 2 ) = MP(Q 2 ) — M£(Q 2 ), together with statistical and systematic uncertainties. 



Q [(GcV/c) 2 ] 


Mj (Q z ) x 10 1 


Ml (Q 2 ) x 10 2 


(Q ) x 10 


(Q z ) x 10 A 


1.025 


6.1 ± 0.6 ± 1.0 


1.74 ± 0.21 ± 0.12 


9.39 ± 0.53 ± 0.12 


7.78 ± 1.01 ± 1.66 


1.075 


6.0 ± 0.6 ± 1.0 


1.67 ± 0.20 ± 0.11 


8.56 ± 0.48 ± 0.11 


6.48 ± 0.84 ± 1.38 


1.125 


6.0 ± 0.6 ± 1.0 


1.62 ± 0.19 ± 0.11 


7.91 ± 0.44 ± 0.11 


5.58 ± 0.73 ± 1.19 


1.175 


5.9 ± 0.6 ± 1.0 


1.57 ± 0.19 ± 0.11 


7.38 ± 0.41 ± 0.11 


4.92 ± 0.64 ± 1.05 


1.225 


5.9 ± 0.6 ± 1.0 


1.52 ± 0.18 ± 0.10 


6.94 ± 0.39 ± 0.10 


4.42 ± 0.57 ± 0.94 


1.275 


5.8 ± 0.6 ± 1.0 


1.48 ± 0.18 ± 0.10 


6.57 ± 0.37 ± 0.10 


4.03 ± 0.52 ± 0.86 


1.325 


5.8 ± 0.6 ± 1.0 


1.45 ± 0.17 ± 0.10 


6.25 ± 0.35 ± 0.10 


3.71 ± 0.48 ± 0.79 


1.375 


5.8 ± 0.6 ± 1.0 


1.42 ± 0.17 ± 0.10 


5.98 ± 0.34 ± 0.10 


3.45 ± 0.45 ± 0.74 


1.425 


5.7 ± 0.6 ± 1.0 


1.39 ± 0.17 ± 0.09 


5.74 ± 0.32 ± 0.09 


3.23 ± 0.42 ± 0.69 


1.475 


5.7 ± 0.6 ± 1.0 


1.36 ± 0.16 ± 0.09 


5.53 ± 0.31 ± 0.09 


3.04 ± 0.40 ± 0.65 


1.525 


5.7 ± 0.6 ± 1.0 


1.34 ± 0.16 ± 0.09 


5.35 ± 0.30 ± 0.09 


2.88 ± 0.37 ± 0.62 


1.575 


5.7 ± 0.6 ± 0.9 


1.31 ± 0.16 ± 0.09 


5.18 ± 0.29 ± 0.09 


2.75 ± 0.36 ± 0.59 


1.625 


5.6 ± 0.6 ± 0.9 


1.29 ± 0.16 ± 0.09 


5.03 ± 0.28 ± 0.09 


2.62 ± 0.34 ± 0.56 


1.675 


5.6 ± 0.6 ± 0.9 


1.27 ± 0.15 ± 0.09 


4.89 ± 0.27 ± 0.09 


2.52 ± 0.33 ± 0.54 


1.725 


5.6 ± 0.6 ± 0.9 


1.26 ± 0.15 ± 0.08 


4.77 ± 0.27 ± 0.08 


2.42 ± 0.31 ± 0.52 


1.775 


5.6 ± 0.6 ± 0.9 


1.24 ± 0.15 ± 0.08 


4.65 ± 0.26 ± 0.08 


2.33 ± 0.30 ± 0.50 


1.825 


5.5 ± 0.6 ± 0.9 


1.22 ± 0.15 ± 0.08 


4.54 ± 0.26 ± 0.08 


2.25 ± 0.29 ± 0.48 


1.875 


5.5 ± 0.6 ± 0.9 


1.21 ± 0.14 ± 0.08 


4.45 ± 0.25 ± 0.08 


2.18 ± 0.28 ± 0.47 


1.925 


5.5 ± 0.6 ± 0.9 


1.19 ± 0.14 ± 0.08 


4.36 ± 0.24 ± 0.08 


2.12 ± 0.28 ± 0.45 


1.975 


5.5 ± 0.6 ± 0.9 


1.18 ± 0.14 ± 0.08 


4.27 ± 0.24 ± 0.08 


2.06 ± 0.27 ± 0.44 


2.025 


5.4 ± 0.6 ± 0.9 


1.17 ± 0.14 ± 0.08 


4.19 ± 0.24 ± 0.08 


2.00 ± 0.26 ± 0.43 


2.075 


5.4 ± 0.6 ± 0.9 


1.16 ± 0.14 ± 0.08 


4.12 ± 0.23 ± 0.08 


1.95 ± 0.25 ± 0.42 


2.125 


5.4 ± 0.6 ± 0.9 


1.15 ± 0.14 ± 0.08 


4.05 ± 0.23 ± 0.08 


1.91 ± 0.25 ± 0.41 


2.175 


5.4 ± 0.6 ± 0.9 


1.13 ± 0.14 ± 0.08 


3.99 ± 0.22 ± 0.08 


1.86 ± 0.24 ± 0.40 


2.225 


5.4 ± 0.6 ± 0.9 


1.12 ± 0.13 ± 0.08 


3.93 ± 0.22 ± 0.08 


1.82 ± 0.24 ± 0.39 


2.275 


5.4 ± 0.6 ± 0.9 


1.11 ± 0.13 ± 0.07 


3.87 ± 0.22 ± 0.07 


1.79 ± 0.23 ± 0.38 


2.325 


5.3 ± 0.6 ± 0.9 


1.11 ± 0.13 ± 0.07 


3.82 ± 0.21 ± 0.07 


1.76 ± 0.23 ± 0.38 


2.375 


5.3 ± 0.6 ± 0.9 


1.10 ± 0.13 ± 0.07 


3.78 ± 0.21 ± 0.07 


1.73 ± 0.23 ± 0.37 


2.425 


5.3 ± 0.6 ± 0.9 


1.09 ± 0.13 ± 0.07 


3.74 ± 0.21 ± 0.07 


1.70 ± 0.22 ± 0.36 


2.475 


5.3 ± 0.6 ± 0.9 


1.08 ± 0.13 ± 0.07 


3.70 ± 0.21 ± 0.07 


1.68 ± 0.22 ± 0.36 


2.525 


5.3 ± 0.6 ± 0.9 


1.08 ± 0.13 ± 0.07 


3.66 ± 0.21 ± 0.07 


1.66 ± 0.22 ± 0.35 


2.575 


5.3 ± 0.6 ± 0.9 


1.07 ± 0.13 ± 0.07 


3.62 ± 0.20 ± 0.07 


1.63 ± 0.21 ± 0.35 


2.625 


5.3 ± 0.6 ± 0.9 


1.06 ± 0.13 ± 0.07 


3.58 ± 0.20 ± 0.07 


1.61 ± 0.21 ± 0.34 


2.675 


5.3 ± 0.6 ± 0.9 


1.06 ± 0.13 ± 0.07 


3.55 ± 0.20 ± 0.07 


1.59 ± 0.21 ± 0.34 


2.725 


5.2 ± 0.6 ± 0.9 


1.05 ± 0.13 ± 0.07 


3.52 ± 0.20 ± 0.07 


1.57 ± 0.20 ± 0.33 


2.775 


5.2 ± 0.6 ± 0.9 


1.04 ± 0.13 ± 0.07 


3.48 ± 0.20 ± 0.07 


1.55 ± 0.20 ± 0.33 


2.825 


5.2 ± 0.6 ± 0.9 


1.04 ± 0.12 ± 0.07 


3.45 ± 0.19 ± 0.07 


1.53 ± 0.20 ± 0.33 


2.875 


5.2 ± 0.6 ± 0.9 


1.03 ± 0.12 ± 0.07 


3.42 ± 0.19 ± 0.07 


1.51 ± 0.20 ± 0.32 


2.925 


5.2 ± 0.6 ± 0.9 


1.03 ± 0.12 ± 0.07 


3.40 ± 0.19 ± 0.07 


1.50 ± 0.19 ± 0.32 


2.975 


5.2 ± 0.6 ± 0.9 


1.02 ± 0.12 ± 0.07 


3.37 ± 0.19 ± 0.07 


1.48 ± 0.19 ± 0.32 


3.025 


5.2 ± 0.5 ± 0.9 


1.02 ± 0.12 ± 0.07 


3.34 ± 0.19 ± 0.07 


1.47 ± 0.19 ± 0.31 


3.075 


5.2 ± 0.5 ± 0.9 


1.01 ± 0.12 ± 0.07 


3.32 ± 0.19 ± 0.07 


1.45 ± 0.19 ± 0.31 


3.125 


5.2 ± 0.5 ± 0.9 


1.01 ± 0.12 ± 0.07 


3.29 ± 0.18 ± 0.07 


1.44 ± 0.19 ± 0.31 


3.175 


5.1 ± 0.5 ± 0.9 


1.00 ± 0.12 ± 0.07 


3.27 ± 0.18 ± 0.07 


1.42 ± 0.18 ± 0.30 


3.225 


5.1 ± 0.5 ± 0.9 


1.00 ± 0.12 ± 0.07 


3.24 ± 0.18 ± 0.07 


1.41 ± 0.18 ± 0.30 


3.275 


5.1 ± 0.5 ± 0.9 


0.99 ± 0.12 ± 0.07 


3.22 ± 0.18 ± 0.07 


1.40 ± 0.18 ± 0.30 


3.325 


5.1 ± 0.5 ± 0.9 


0.99 ± 0.12 ± 0.07 


3.20 ± 0.18 ± 0.07 


1.38 ± 0.18 ± 0.29 


3.375 


5.1 ± 0.5 ± 0.8 


0.98 ± 0.12 ± 0.07 


3.18 ± 0.18 ± 0.07 


1.37 ± 0.18 ± 0.29 


3.425 


5.1 ± 0.5 ± 0.8 


0.98 ± 0.12 ± 0.07 


3.16 ± 0.18 ± 0.07 


1.36 ± 0.18 ± 0.29 


3.475 


5.1 ± 0.5 ± 0.8 


0.98 ± 0.12 ± 0.07 


3.14 ± 0.18 ± 0.07 


1.35 ± 0.18 ± 0.29 


3.525 


5.1 ± 0.5 ± 0.8 


0.97 ± 0.12 ± 0.07 


3.12 ± 0.17 ± 0.07 


1.34 ± 0.17 ± 0.28 


3.575 


5.1 ± 0.5 ± 0.8 


0.97 ± 0.12 ± 0.07 


3.10 ± 0.17 ± 0.07 


1.32 ± 0.17 ± 0.28 


3.625 


5.1 ± 0.5 ± 0.8 


0.96 ± 0.12 ± 0.06 


3.08 ± 0.17 ± 0.06 


1.31 ± 0.17 ± 0.28 


3.675 


5.1 ± 0.5 ± 0.8 


0.96 ± 0.11 ± 0.06 


3.06 ± 0.17 ± 0.06 


1.30 ± 0.17 ± 0.28 


3.725 


5.1 ± 0.5 ± 0.8 


0.96 ± 0.11 ± 0.06 


3.04 ± 0.17 ± 0.06 


1.29 ± 0.17 ± 0.28 


3.775 


5.0 ± 0.5 ± 0.8 


0.95 ± 0.11 ± 0.06 


3.03 ± 0.17 ± 0.06 


1.28 ± 0.17 ± 0.27 


3.825 


5.0 ± 0.5 ± 0.8 


0.95 ± 0.11 ± 0.06 


3.01 ± 0.17 ± 0.06 


1.28 ± 0.17 ± 0.27 
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Q 2 [(GcV/c) 2 l 


M|""(Q 2 ) x 10- 2 


M^ n (Q 2 ) x 10~ 2 


MS~ n (Q 2 ) x 10" 3 


Mg- n (Q 2 ) x 10- 3 


3.875 


5.0 ± 0.5 ± 0.8 


0.95 ± 0.11 ± 0.06 


2.99 ± 0.17 ± 0.06 


1.27 ± 0.16 ± 0.27 


3.925 


5.0 ± 0.5 ± 0.8 


0.94 ± 0.11 ± 0.06 


2.98 ± 0.17 ± 0.06 


1.26 ± 0.16 ± 0.27 


3.975 


5.0 ± 0.5 ± 0.8 


0.94 ± 0.11 ± 0.06 


2.96 ± 0.17 ± 0.06 


1.25 ± 0.16 ± 0.27 


4.025 


5.0 ± 0.5 ± 0.8 


0.94 ± 0.11 ± 0.06 


2.95 ± 0.17 ± 0.06 


1.24 ± 0.16 ± 0.26 


4.075 


5.0 ± 0.5 ± 0.8 


0.93 ± 0.11 ± 0.06 


2.93 ± 0.16 ± 0.06 


1.23 ± 0.16 ± 0.26 


4.125 


5.0 ± 0.5 ± 0.8 


0.93 ± 0.11 ± 0.06 


2.92 ± 0.16 ± 0.06 


1.22 ± 0.16 ± 0.26 


4.175 


5.0 ± 0.5 ± 0.8 


0.93 ± 0.11 ± 0.06 


2.90 ± 0.16 ± 0.06 


1.22 ± 0.16 ± 0.26 


4.225 


5.0 ± 0.5 ± 0.8 


0.92 ± 0.11 ± 0.06 


2.89 ± 0.16 ± 0.06 


1.21 ± 0.16 ± 0.26 


4.275 


5.0 ± 0.5 ± 0.8 


0.92 ± 0.11 ± 0.06 


2.87 ± 0.16 ± 0.06 


1.20 ± 0.16 ± 0.26 


4.325 


5.0 ± 0.5 ± 0.8 


0.92 ± 0.11 ± 0.06 


2.86 ± 0.16 ± 0.06 


1.19 ± 0.16 ± 0.25 


4.375 


5.0 ± 0.5 ± 0.8 


0.92 ± 0.11 ± 0.06 


2.85 ± 0.16 ± 0.06 


1.19 ± 0.15 ± 0.25 


4.425 


5.0 ± 0.5 ± 0.8 


0.91 ± 0.11 ± 0.06 


2.84 ± 0.16 ± 0.06 


1.18 ± 0.15 ± 0.25 


4.475 


5.0 ± 0.5 ± 0.8 


0.91 ± 0.11 ± 0.06 


2.82 ± 0.16 ± 0.06 


1.17 ± 0.15 ± 0.25 


4.525 


4.9 ± 0.5 ± 0.8 


0.91 ± 0.11 ± 0.06 


2.81 ± 0.16 ± 0.06 


1.17 ± 0.15 ± 0.25 


4.575 


4.9 ± 0.5 ± 0.8 


0.91 ± 0.11 ± 0.06 


2.80 ± 0.16 ± 0.06 


1.16 ± 0.15 ± 0.25 


4.625 


4.9 ± 0.5 ± 0.8 


0.90 ± 0.11 ± 0.06 


2.79 ± 0.16 ± 0.06 


1.15 ± 0.15 ± 0.25 


4.675 


4.9 ± 0.5 ± 0.8 


0.90 ± 0.11 ± 0.06 


2.78 ± 0.16 ± 0.06 


1.15 ± 0.15 ± 0.25 


4.725 


4.9 ± 0.5 ± 0.8 


0.90 ± 0.11 ± 0.06 


2.76 ± 0.16 ± 0.06 


1.14 ± 0.15 ± 0.24 


4.775 


4.9 ± 0.5 ± 0.8 


0.90 ± 0.11 ± 0.06 


2.75 ± 0.15 ± 0.06 


1.14 ± 0.15 ± 0.24 


4.825 


4.9 ± 0.5 ± 0.8 


0.89 ± 0.11 ± 0.06 


2.74 ± 0.15 ± 0.06 


1.13 ± 0.15 ± 0.24 


4.875 


4.9 ± 0.5 ± 0.8 


0.89 ± 0.11 ± 0.06 


2.73 ± 0.15 ± 0.06 


1.13 ± 0.15 ± 0.24 


4.925 


4.9 ± 0.5 ± 0.8 


0.89 ± 0.11 ± 0.06 


2.72 ± 0.15 ± 0.06 


1.12 ± 0.15 ± 0.24 


4.975 


4.9 ± 0.5 ± 0.8 


0.89 ± 0.11 ± 0.06 


2.71 ± 0.15 ± 0.06 


1.11 ± 0.14 ± 0.24 


5.025 


4.9 ± 0.5 ± 0.8 


0.88 ± 0.11 ± 0.06 


2.70 ± 0.15 ± 0.06 


1.11 ± 0.14 ± 0.24 


5.075 


4.9 ± 0.5 ± 0.8 


0.88 ± 0.11 ± 0.06 


2.69 ± 0.15 ± 0.06 


1.10 ± 0.14 ± 0.24 


5.125 


4.9 ± 0.5 ± 0.8 


0.88 ± 0.11 ± 0.06 


2.68 ± 0.15 ± 0.06 


1.10 ± 0.14 ± 0.23 


5.275 


4.9 ± 0.5 ± 0.8 


0.87 ± 0.10 ± 0.06 


2.65 ± 0.15 ± 0.06 


1.08 ± 0.14 ± 0.23 


5.325 


4.9 ± 0.5 ± 0.8 


0.87 ± 0.10 ± 0.06 


2.64 ± 0.15 ± 0.06 


1.08 ± 0.14 ± 0.23 


5.375 


4.9 ± 0.5 ± 0.8 


0.87 ± 0.10 ± 0.06 


2.64 ± 0.15 ± 0.06 


1.07 ± 0.14 ± 0.23 


5.475 


4.9 ± 0.5 ± 0.8 


0.87 ± 0.10 ± 0.06 


2.62 ± 0.15 ± 0.06 


1.07 ± 0.14 ± 0.23 


5.525 


4.8 ± 0.5 ± 0.8 


0.86 ± 0.10 ± 0.06 


2.61 ± 0.15 ± 0.06 


1.06 ± 0.14 ± 0.23 


5.625 


4.8 ± 0.5 ± 0.8 


0.86 ± 0.10 ± 0.06 


2.59 ± 0.15 ± 0.06 


1.05 ± 0.14 ± 0.22 


5.675 


4.8 ± 0.5 ± 0.8 


0.86 ± 0.10 ± 0.06 


2.59 ± 0.15 ± 0.06 


1.05 ± 0.14 ± 0.22 


5.725 


4.8 ± 0.5 ± 0.8 


0.86 ± 0.10 ± 0.06 


2.58 ± 0.14 ± 0.06 


1.04 ± 0.14 ± 0.22 


5.955 


4.8 ± 0.5 ± 0.8 


0.85 ± 0.10 ± 0.06 


2.54 ± 0.14 ± 0.06 


1.03 ± 0.13 ± 0.22 


6.915 


4.7 ± 0.5 ± 0.8 


0.82 ± 0.10 ± 0.06 


2.42 ± 0.14 ± 0.06 


0.96 ± 0.12 ± 0.21 


7.267 


4.7 ± 0.5 ± 0.8 


0.81 ± 0.10 ± 0.05 


2.38 ± 0.13 ± 0.05 


0.94 ± 0.12 ± 0.20 


7.630 


4.7 ± 0.5 ± 0.8 


0.80 ± 0.10 ± 0.05 


2.34 ± 0.13 ± 0.05 


0.92 ± 0.12 ± 0.20 


8.021 


4.7 ± 0.5 ± 0.8 


0.79 ± 0.10 ± 0.05 


2.31 ± 0.13 ± 0.05 


0.91 ± 0.12 ± 0.19 


8.847 


4.6 ± 0.5 ± 0.8 


0.78 ± 0.09 ± 0.05 


2.24 ± 0.13 ± 0.05 


0.87 ± 0.11 ± 0.19 


9.775 


4.6 ± 0.5 ± 0.8 


0.76 ± 0.09 ± 0.05 


2.17 ± 0.12 ± 0.05 


0.84 ± 0.11 ± 0.18 


10.267 


4.6 ± 0.5 ± 0.8 


0.76 ± 0.09 ± 0.05 


2.14 ± 0.12 ± 0.05 


0.83 ± 0.11 ± 0.18 


10.762 


4.6 ± 0.5 ± 0.8 


0.75 ± 0.09 ± 0.05 


2.12 ± 0.12 ± 0.05 


0.81 ± 0.11 ± 0.17 


11.344 


4.5 ± 0.5 ± 0.8 


0.74 ± 0.09 ± 0.05 


2.08 ± 0.12 ± 0.05 


0.80 ± 0.10 ± 0.17 


12.580 


4.5 ± 0.5 ± 0.7 


0.73 ± 0.09 ± 0.05 


2.03 ± 0.11 ± 0.05 


0.77 ± 0.10 ± 0.16 


13.238 


4.5 ± 0.5 ± 0.7 


0.72 ± 0.09 ± 0.05 


2.00 ± 0.11 ± 0.05 


0.76 ± 0.10 ± 0.16 


14.689 


4.4 ± 0.5 ± 0.7 


0.71 ± 0.08 ± 0.05 


1.95 ± 0.11 ± 0.05 


0.73 ± 0.09 ± 0.16 


17.108 


4.4 ± 0.5 ± 0.7 


0.69 ± 0.08 ± 0.05 


1.87 ± 0.11 ± 0.05 


0.70 ± 0.09 ± 0.15 


19.072 


4.3 ± 0.5 ± 0.7 


0.67 ± 0.08 ± 0.05 


1.82 ± 0.10 ± 0.05 


0.67 ± 0.09 ± 0.14 


20.108 


4.3 ± 0.5 ± 0.7 


0.67 ± 0.08 ± 0.04 


1.80 ± 0.10 ± 0.04 


0.66 ± 0.09 ± 0.14 


21.097 


4.3 ± 0.5 ± 0.7 


0.66 ± 0.08 ± 0.04 


1.78 ± 0.10 ± 0.04 


0.66 ± 0.09 ± 0.14 


24.259 


4.3 ± 0.5 ± 0.7 


0.65 ± 0.08 ± 0.04 


1.73 ± 0.10 ± 0.04 


0.63 ± 0.08 ± 0.13 


26.680 


4.2 ± 0.5 ± 0.7 


0.64 ± 0.08 ± 0.04 


1.69 ± 0.09 ± 0.04 


0.61 ± 0.08 ± 0.13 


32.500 


4.2 ± 0.4 ± 0.7 


0.62 ± 0.07 ± 0.04 


1.62 ± 0.09 ± 0.04 


0.58 ± 0.08 ± 0.12 


34.932 


4.2 ± 0.4 ± 0.7 


0.61 ± 0.07 ± 0.04 


1.60 ± 0.09 ± 0.04 


0.57 ± 0.07 ± 0.12 


36.750 


4.1 ± 0.4 ± 0.7 


0.61 ± 0.07 ± 0.04 


1.58 ± 0.09 ± 0.04 


0.57 ± 0.07 ± 0.12 


43.970 


4.1 ± 0.4 ± 0.7 


0.60 ± 0.07 ± 0.04 


1.53 ± 0.09 ± 0.04 


0.54 ± 0.07 ± 0.12 


47.440 


4.1 ± 0.4 ± 0.7 


0.59 ± 0.07 ± 0.04 


1.51 ± 0.08 ± 0.04 


0.54 ± 0.07 ± 0.11 


64.270 


4.0 ± 0.4 ± 0.7 


0.57 ± 0.07 ± 0.04 


1.43 ± 0.08 ± 0.04 


0.50 ± 0.06 ± 0.11 


75.000 


4.0 ± 0.4 ± 0.7 


0.55 ± 0.07 ± 0.04 


1.39 ± 0.08 ± 0.04 


0.48 ± 0.06 ± 0.10 


86.000 


3.9 ± 0.4 ± 0.7 


0.55 ± 0.07 ± 0.04 


1.36 ± 0.08 ± 0.04 


0.47 ± 0.06 ± 0.10 


97.690 


3.9 ± 0.4 ± 0.7 


0.54 ± 0.06 ± 0.04 


1.33 ± 0.07 ± 0.04 


0.46 ± 0.06 ± 0.10 
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Figure 10: Ratio of extracted neutron to proton leading twist moments as a function ofn (solid), 
compared with the structure function ratio (dashed) evaluated using the PDFs of Ref. ]3Ctf . as 
a function of (x) n (top axis). The horizontal lines denote the asymptotic x — > 1 values for the 
structure function moments of 1/4 \34\ an d 3/7 f.Vffi . 



in the quark masses up to the physical point, while in Ref. [HEj the extrapolation includes 
the effects of meson loops and intermediate A(1232) resonance which are known to 
provide important non-analytic terms present in chiral perturbation theory. 

Using Eq. (|38j) the results of Refs. |37| I38| I39| are translated in terms of nonsinglet 
p — n moments, which are reported in Table |I] together with our corresponding results 
from Table El It can be seen that the results of Refs. [S3 EH]) which adopt a naive linear 
extrapolation in the quark masses up to the chiral point, are larger than our findings 
both at n = 2 and n = 4. Note that, by combining all the uncertainties, the deviations 
between our data and the lattice results of Refs. j33 OH] are » 3.3a at n = 2, but only 
~ 1.6a at n = 4. The results of Ref. [22], which includes the effects of meson loops and 
intermediate A(1232) resonance in the chiral extrapolation, are in excellent agreement 
with our numbers for both n = 2 and n = 4. It should be mentioned that the role of 
finite volume effects at small quark masses still remains to be investigated. 

Table 4: Comparison of the nonsinglet leading twist moments of the F2 structure function given 
in Tabled with the lattice predictions from Refs. H%1 HQ/ at the scale Q 2 = 4 GeV 2 . 



n 


This work 


Ref. 37 


Ref. 38 


Ref. 39 


2 


0.050 (5) (8) 


0.091 (8) 


0.082 (3) 


0.059 (8) 


4 


0.0094 (11) (6) 


0.030 (16) 


0.023 (7) 


0.009 (3) 
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Figure 11: Ratio of the extracted neutron to proton leading twist moments (triangles) as a 
function of n, compared with the corresponding PDF parameterization of Ref. (solid). The 
hatched area represents our systematic errors. The dashed and dotted curves correspond to the 
moment ratio evaluated after modifying the d-quark distribution according to Eq. |ff7| ) with m = 1 
and m — 4 7 respectively. 

4 Isospin dependence of the Higher Twist contribution 

Once the leading twist is settled and checked, one can study the isospin dependence of the 
higher twist (HT) contribution. Previous studies were performed in x-space in Refs. |4()l 
I41j . These analyses showed that in the region x < 0.7 HT terms are approximately the 
same for the proton and the deuteron, while at larger x some isospin dependence was 
seen in Ref. |41j . but not in Ref. |4L)j (although the experimental errors in Ref. |40| were 
much larger than the difference between the HTs in the proton and in the deuteron). 

In the present paper we consider the total HT contribution in moment space, simply 
defined as the difference between the total and leading twist moments of F2, 

HT^(Q 2 ) = M^ tot) (Q 2 ) - M^(Q 2 ) , (39) 

where the leading twist moment M^(Q 2 ) is defined as in Eq. (J2J). 

Let us briefly recall that our approach is characterized by some relevant features 
which are essential for a reliable extraction of the HT contribution. In particular: (i) 
the target mass corrections are removed from the total moments through the use of 
the Nachtmann definition; (ii) effective HT anomalous dimensions are extracted phe- 
nomenologically; (iii) the effects of Sudakov logarithms are taken into account through 
Soft Gluon Resummation; and (iv) the running coupling constant a s (M|) is consis- 
tently extracted from the proton moments at high Q 2 |42j . giving the value a s (M|) = 




2 4 6 8 1( 

n 
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0.1188 ± 0.0010 s tat i 0.0014 syst , which is consistent with the latest world average value 
a s (M|) = 0.1187± 0.0020 As a result the total higher twist contribution HT^ (Q 2 ) 
is extracted reliably and with good precision both for the proton and deuteron Ej for 
Q 2 >lGeV 2 . 

In Fig. ^]the total HT term in the proton moments, HT^(Q 2 ), is compared with 
that in the deuteron, HT^ (Q 2 ), after applying to the latter the same nuclear corrections 
used for the leading twist, namely of Eq. (|14|). Within the quoted uncertainties the 
higher twists appear to be the same in the proton and deuteron (per nucleon) moments. 
Assuming that possible additional nuclear effects are small with respect to the moment 
uncertainties (at least for Q 2 > 1 GeV 2 ), the total HT contribution appears to be isospin- 
independent. This means that for the isovector combination p — n our total HT term is 
consistent with zero within the uncertainties. 

Based on the above results we now speculate about two possible scenarios in which 
our findings can be interpreted: 

1. isospin-independent q — q correlations. 

2. dominance of u — d correlations with respect to u — u and d — d ones. 

Assuming that HTs are dominated by coupling of the virtual boson with correlated qq 
pairs, the neutron to proton ratio of HTs is expected to be 

HT n = 2(ud)\ + (dd)l 

HTP 2{ud)l + (uu)f { ' 

If q — q correlations are isospin-independent, the HT ratio (|40[) becomes equal to 1/3, a 
value which is not supported by the present analysis. 

The second scenario can be better understood in terms of a diquark model (see 
|44j and references therein), which represent a simple, natural way to account for the 
so-called lack of missing nucleon resonances. The correlation between two quarks in a 
+ color antitriplet state is expected to be of particular importance jSlEHl- Indeed, 
such a configuration is very likely in the nucleon because the correlation between two 
quarks in a color antitriplet state forming a scalar, isoscalar diquark is as strong as the 
corresponding qq correlation in the pion, as is known from both QCD sum rules |4fij 
and instanton calculations [IT]- The most dramatic effect of scalar, isoscalar diquark 
formation occurs in strangeness changing weak decays at low energies, leading to the 
explanation of the huge AI = 1/2 enhancement present in these processes |48j . 

In the quark-diquark picture of the nucleon diquarks are expected to be spatially 
extended objects and therefore it is natural to consider that HT terms in lepton-nucleon 
DIS originate mainly from the elastic scattering of the virtual boson off diquarks, rep- 
resenting correlated pairs of quarks. Many applications of the quark-diquark model to 
nucleon observables indicate that scalar, isoscalar diquarks are more abundant and have 
a size smaller than the one of vector, isovector diquarks 44 . This means that u — d 
correlations are expected to dominate over u — u and d — d ones, and therefore the HT 
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Figure 12: TotaZ higher twist contribution to the moments with n = 2, 4, 6, 8 together with their 
statistical (bars) and systematic (left- and across-hatched areas) uncertainties: triangles show the 
higher twist contribution in the proton moments and the stars represent the same contribution 
in the deuteron moments corrected by the Fermi motion factor given in Eq. The left- 

and cross-hatched areas refer to the case of proton and deuteron, respectively. 



ratio (|4Uj) becomes equal to 1, a value which is fully consistent with our findings. A clear- 
cut determination of the spin of the dominant u — d diquark requires the simultaneous 
investigation of the HTs in the longitudinal channel. 

For completeness we should mention that the dominance of isoscalar over isovector 
higher-twist contributions in unpolarized structure functions is qualitatively predicted by 
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the 1/N C expansion of QCD, as was recently pointed out in the context of a calculation 
of twist-4 matrix elements in the chiral quark-soliton model |49| . 

Before closing this Section, we want to point out that, if HTs are isospin-independent, 
then the nuclear correction can be obtained directly from the ratio between higher twist 
contributions in the deuteron and proton moments. In Fig. ^3 such a ratio is compared 
with the nuclear correction factor given by Eq. (|14jl. The agreement is reasonable within 
the statistical and systematic error bars. 
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Figure 13: Ratio of the deuteron to the proton higher twist contribution at Q 2 — 3 GeV 2 with its 
statistical and systematic uncertainties represented by bars and hatched area, respectively. The 
solid curve represents the nuclear correction factor defined in Eq. 



5 Conclusions and future directions 

We have extracted the leading twist contribution to first few moments (n = 2 — 10) of 
the neutron structure function Fi in the Q 2 range from 1 up to 100 GeV 2 , by combining 
proton and deuteron measurements over a huge range of kinematics. We have paid 
particular attention to the issue of nuclear effects in the deuteron, which are found to 
be increasingly important for the higher moments. The model dependence arising from 
our incomplete knowledge of the nuclear corrections, including the high-momentum tail 
of the deuteron wave function, as well as nucleon off-shell effects, were estimated and 
included in the systematic errors. From our analysis the following conclusions can be 
made: 

• the extracted moments of the neutron structure function F2 are consistent at large 
Q 2 with the predictions of several PDFs available in the literature; 
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• the ratio of neutron to proton structure functions is determined from the ratio of 
the corresponding moments up to x — 0.70; we find FQ jF^ [x = 0.70) = 0.34 ± 
0.12 sta t i 0.13 sys t. Our results are consistent with the asymptotic limit 1/4 at 
x = 1, which originates from the dominance of soft, non-perturbative physics at 
large x. Nevertheless the alternative limiting value of 3/7, derived from helicity 
conservation arguments, is not excluded at all, but the possible enhancement of 
the d-quark distribution should be confined to the region x > 0.7; 

• the nonsinglet, isovector p — n combination of the leading twist moments is con- 
structed and compared with available lattice results. We have found that our data 
are larger than the lattice values of Refs. |37| I38j . extrapolated linearly in the quark 
masses down to the physical point, by ~ 3. 3d for the second moment, but only by 
~ 1.6(7 in case of the fourth moment. Our results are in excellent agreement with 
those of Ref . 23 > where the chiral extrapolation is performed taking into account 
the effects of non-analytic terms due to meson loops and intermediate A(1232) 
resonance; 

• the total contribution of higher twists is found to be isospin independent, which 
implies that in the isovector combination p — n the higher twists are consistent 
with zero within the uncertainties. The isospin independence is expected in a 
quark-diquark picture of the nucleon assuming the dominance of isoscalar u — d 
diquarks. 

The determination of the large- £ behavior of the d-quark distribution in the nucleon 
is still one of the main unsolved problems in partonic physics. The use of deuteron 
inclusive DIS data to extract the neutron to proton structure function ratio produces 
unavoidably a model dependence. The strategy used in this paper has its own practical 
limitations which are due to the increasing uncertainties of the nuclear corrections as 
the order of the moments considered increases. It appears that moments of order n = 12 
may be still accessible by improving measurements of the structure function Fi at large 
x. As we have shown, at n = 12 one can probe values of the Bjorken-x as large as 0.75, 
which represents therefore an upper limit. 

Future directions for the experimental determination of the (i-quark distribution at 
high x have been recently discussed in the literature [SOI 03 EH 023 • One of the strategies, 
which is currently under experimental investigation at Jefferson Lab is the use of 
the tagged semi-inclusive process off the deuteron (SUJ03, which requires the detection 
of a low momentum proton (p < 150MeV/c). An alternative method which has been 
discussed, tries to exploit the mirror symmetry of A = 3 nuclei to reduce nuclear effects 
in the extraction of the FQ / F^ ratio |521 153j . 

In the semi-inclusive deuteron method one can tag the momentum of the struck 
neutron by detecting the slow recoiling proton; in this way it is possible to select initial 
deuteron configurations in which the two nucleons are far apart, so that the struck 
nucleon can be considered as free and the uncertainties due to the low-momentum part 
of the deuteron wave function are minimal |50| I51j . The neutron structure function can 
then be extracted directly from the semi-inclusive deuteron cross section. 
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We want to point out that an important improvement of the tagged strategy may be 
achieved by adding the detection of slow recoiling neutrons. In this way the n/p structure 
function ratio can be directly obtained from the ratio of the corresponding semi-inclusive 
cross sections. Most of the nuclear corrections, e.g. Fermi motion, off-shellness and FSIs, 
are expected to largely cancel out in the ratio |5U1 15 lj . Thus, from the measurements 
of both semi-inclusive processes D(e,e'p)X and D{e,e'n)X the shape of the d-quark 
distribution may be experimentally investigated up to and hopefully beyond x « 0.75. 
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